JPRS-UEN-91-92 
8 FEBRUARY 1991 











JPRS Report— 








Science & 
Technology 





USSR: Engineering & 
Equipment 











JPRS-UEQ-91-102 
8 FEBRUARY 1991 


SCIENCE & TECHNOLOGY 
USSR: ENGINEERING & EQUIPMENT 


CONTENTS 


MECHANICS OF GASES, LIQUIDS, SOLIDS 


Magnetoelastic Surface Waves in Half-Space With Surface 
Electric Current 
[M. V. Belubekyan; IZVESTIYA AKADEMII NAUK 
ARMYANSKOY SSR: SERIYA MEKHANIKA, Vol 43 No 1, 1990]...... 


Diffraction of Shock Wave Around Wall Making Obtuse Angle 
[L. D. Azatyan, A. G. Bagdoyev; IZVESTIYA AKADEMII NAUK 
ARMYANSKOY SSR: SERIYA MEKHANIKA, Vol 43 No 1, 1990]...... 


Frequencies of Free Oscillations of Transversely Loaded Plates 
[V. Ts. Gnuni; IZVESTIYA AKADEMII NAUK ARMYANSKOY SSR: 
SERIYA MEKHANIKA, Vol 43 No 1, 1990]... cc ccc ccc cecceveee 


Design of Optimum Circular Plate of Maximum Stiffness Made of 
Composite Material 
[S. A. Grigoryan; IZVESTIYA AKADEMII NAUK ARMYANSKOY SSR: 
SERIYA MEKHANIKA, Vol 43 No 1, 1990]... cc cece cece cree eens 


Wave Propagation in Radially Inhomogeneous Elastic Ball 
[S. G. Saakyan; IZVESTIYA AKADEMII NAUK ARMYANSKOY SSR: 
SERIYA MEKHANIKA, Vol 43 No 1, 1990]... .. ccc erence enc cccens 








Introduction of Rigid Cylinder Into Cylinder Made of Material 
Conforming to Power Hardening Law 
(Zh. G. Apikyan; IZVESTIYA AKADEMII NAUK ARMYANSKOY SSR: 
SERIYA MEKHANTIKA, Vol 43 No 1, 1990] 


: 60:0 ©% €054566004 45400084 34 
Electrical Transient Processes in the Vapor-Gas Shell Under 
Conditions of Electrolytic Anodic Heating With Rectangular 
Voltage Pulses 
[S. I. Galanin, V. I. Ganchar, et al.; ELEKTRONNAYA 
OBRABOTKA MATERIALOV, No 2, 1990]... ccccccccceccecccuvuces 43 
Investigation of an Underwater Electric Explosion Initiated 
by Ferromagnetic Particles. I. Physical Modeling 
[L. Z. Boguslavskiy, L. M. Bondarets, et al.; 
ELEKTRONNAYA OBRABOTKA MATERTALOV, No 2, 1990]..........06. 48 





MECHANICS OF GASES, LIQUIDS, SOLIDS 


UDC 539.3:534.21 
Magnetoelastic Surface Waves in Half-Space With Surface Electric Current 


917FO0082A IZVESTIYA AKADEMII NAUK ARMYANSKOY SSR: SERIYA MEKHANIKA in 
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[Article by M. V. Belubekyan, Institute of Mechanics, ArmSSR Academy of 
Sciences ] 


[Text] The earliest research on magnetoelastic surface waves in the 
presence of a constant magnetic field was done by S. Kaliskiy |1]. This 
work has subsequently been continued by many researchers. 


In contrast to other papers, this one considers the case where a magnet- 
ic field is set up by electric current flowing over the surface of a 
half-space. Studies of the propagation of elastic waves in the presence 
of surface current may have application in nondestructive testing of 
materials [2]. 


1. Let an elastic half-space occupy region x, 2 QO in rectangular coor- 
dinate system (x,, Xz, X3). Electric current with surface density J, 
flows along boundary x, = O in the direction of axis Ox,. In the ini- 
tial (unperturbed) state, the magnetic fields caused by J, will be 
constant in regions x, > O and x, < 0, undergoing a ¢ciscontinuity (jump) 
on surface x, = 0. The magnetic field intensities must be defined to 
satisfy boundary condition 


—s 


_ — N  4z , f 
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Here HS°?, Ho are magnetic field intensity vectors where x, < 0, x, > 0 
respectively, 1, j, k are unit vectors in directions x,, x, and x, 
respectively. 


The following three versions are possible for definition of a magnetic 
field satisfying these conditions: 
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A model for versions (1.2) and (1.3) is a medium in the form of a layer 
of rather considerable thickness with current flowing over its surfaces 


with surface density J,. Where currents over the surfaces have the same 
direction, the magnetic field cannot penetrate into the layer, and we 
arrive at version (1.2). Version (1.3) is the result in the case where 


surface currents are oppositely directed. 


According to (1.1)-(1.3, the components of the Maxwell tensors for the 
corresponding versions take the form (only non-zero components are 
given): 
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Here and below, the permeabilities of the material of the half-space 
x, 2 9 and of the medium outside of the half-space are assumed to be 
equal to unity. 


In the first version (1.4), the components of the Maxwell tensor are 
continuous on boundary x. = O, and therefore the initial stressed state 
is zero (of, = 0). In the second (1.5) and third (1.6) versions, compo- 
nent T,3:,. undergoes a discontinuity on boundary x, = 0, giving rise to 
initial stresses. 


To determine the initial stressed state, it is necessary to satisfy 
boundary conditions (x. = 0) 
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The stressed state also depends on conditions as x, > +, x, > +, The 
simplest stressed state results when conditions 


lime},=9, lim? = U (1.8) 
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are satisfied. 


In this context, all components of the stress tensor are equal to zero 
except of., which is defined as follows 
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(for the second and third versions respectively). 


Formulation of the problem of wave propagation in the given three ver- 
Sions of the initial state will depend on the physical properties of the 
material of the half-space. Here problems are possible where the ma- 
terial is either a dielectric, or an ideal conductor, or a superconduc- 
tor. 


2. Of considerable significance in derivation of boundary conditions is 
allowance for the behavior of the normal to the perturbed surface of the 
half-space [3]. 


Let the perturbed surface of the half-space be defined by the equation 
Ny Uy(X,. O, X,. 0) (2.1) 


Then the approximate (linearized) expression for the normal to the 
perturbed surface takes the form 


te Oat, Oe (2.2) 
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To derive linearized boundary conditions on surface x, = 0, we will 
start from the following general conditions; when x. = 0 
= : \ | Ou 7 
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Using expressions 
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for the perturbed state, and linearizing condition (2.3) with allowance 
for smallness of perturbations and with consideration of (1.1-1.6), 
(1.9), (2.2), we get the boundary conditions on the surface of the 
half-space. 


In all three versions,. boundary conditions for components of the per- 
turbed electromagnetic field are the same, and take the form 
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Conditions for stresses for versions (1.1)-(1.3) in the corresponding 
sequence are obtained in the form 
a f H, H H 
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In (2.7)-(2.9) and below, the primes of o}, have been omitted. 


Let us take a look at some special problems 
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Let us consider the plane problem where displacements and components of 


the perturbed electromagnetic field do not 


depend on coordinate x,,. 





We give the equations and boundary conditions of the problem. The 
equations of the theory of elasticity in displacements 
Pit, a1, arr, ru, 
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According to (2.6) and (2.7), 
take the form 


the required 


boundary conditions at x, = 0 
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Representing the sought functions in the form 
I(X,. XC) file,) exp ((ot--R.x,) 


and following the scheme of solution of the problem of Rayleigh surface 
waves [4], we get an equation that determines the velocity of a surface 


wave (it is assumed that ¢€ = 1, w*/c* « 1) 
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For x, = 0, equation (3.5) coincides with the Rayleigh equation. Just 


as for the Rayleigh equation, it is shown that equation (3.5) has only 
one real root satisfying condition 0 < n < 1. 


Let us give the equations of the corresponding antiplanar problem. 


The equation of the theory of elasticity in region x, > O takes the form 
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The equations of electrodynamics in region x, > 0 
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Boundary conditions at x, = O take the form 
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We have the task of finding the solutions of equations (3.6)-(3.8) that 
satisfy boundary conditions (3.9) and conditions of damping as x. > +. 








In solving the problem in this formulation, we come to the necessity of 
an additional boundary condition at x, = 0. This condition is given in 
the form 
4- Ou 
iG — 3 n 
A —f = — J, —4+ (3.10) 
c OX, 


which is due to the jump in the tangential component of the magnetic 
Field in consequence of a change in the direction of the surface cur- 
rent. 


With allowance for additional condition (3.10), the problem is complete- 
ly defined. Solution of the problem shows that there is a purely shear 
surface wave with velocity defined by the formula 


H- \3 
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4. Let the initial magnetic field be assigned in form (1.3) (third 
version). It is assumed that the ideal conductor model is applicable to 
the material of medium x, > 0. It is further assumed that the sought 
values do not depend on coordinate x,. In this case, the problem of 
plane deformation leads to solution of the following equations (V2 = 


Ho/tp): 
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The solution of equation (4.1) with condition (4.2) and the condition of 
damping at infinity (x, > ~) shows that the condition of existence of a 
surface wave is the same as in the Rayleigh problem (0 < n < 1). The 
equation that defines the velocity of the surface wave takes the form 


(2— 3,2—7,) {(1—68) 12 — 3(—1/ 2 46 - 5,32) -- (1-9, 3) 4h = 
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{Investigation of equation (4.3) relative to the velocity of the surface 
wave n shows that the equation has two real roots in interval (0, 1). 
One root corresponds to the velocity of the ordinary Rayleigh wave with 
allowance for the influence of the magnetic field. The other root shows 
the presence of a surface wave with a velocity much lower than that of 
the Rayleigh wave. This velocity can be computed from the following 
approximate formula: 





2z(/ +(j;) 
The problem outlined in this section has also been considered in [5-7]. 
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{Text] The paper considers the problem of diffraction of a strong shock 
wave around an elastic wall that forms an oblique angle. Lighthill used 
Chaplygin transformation to study the analogous problem and get an 
analytical solution in the case of a solid wall [2]. 


§1. Formulation of the Boundary Value Problem 
Let a shock wave front move at velocity U, through a quiescent gas and 


impact at time t # O on a corner with sides that form a small angle e 
with a break. Fig. 1 shows the pattern of motion for some t # 0. The 




















lo: 


parameters of homogeneous flow (1) behind the incident shock wave are 
determined from conditions of ist 













(designation of parameters is conventional). 


Let us determine the flow in region (2) behind the reflected shock wave. 
In the case of an elastic wall, assuming that motion in the wall is 
homogeneous, we get in the elastic region 


where A, uw are the elastic Lamé constants, p., is the density of the 
elastic body, w is normal displacement of the elastic medium, and a is 
the velocity of propagation of the elastic longitudinal wave. On bound- 
ary x = O from the condition of concomitant motion of the wall and 


nearby particles of fluid (gas) we nave 








where n is the normal to the wall, and the solution in the elastic 
region is w = f(x + at). Writing the conditions of continuity 
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we can get two conditions on the wall 
Qrn=af (at), %a*t '(ath=—P, 
whence we have 
fo Pi 
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where 8 = 9,3c./Pcra- 


Thus, accounting for a one-dimensional model in an elastic medium gives 
a boundary condition in the form of a Winkler base. Flow parameters in 
region (2) are determined from relations for an oblique shock wave [1]. 
Eliminating P,, p,» from these expressions, with allowance for (1.2), we 
find a cubic equation in order e« for determining the normal velocity V, 
of the reflected shock 
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Then in the same approximation from expressions for an oblique shock 
wave we get the values of the remaining parameters for a plane reflected 
shock 
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where e€' is the angle of reflection determined from the condition of 
inseparability of the incident and reflected shock waves at point M 


at (1.5) 





In the region of nonuniform flow (3), we will give all quantities the 
subscript "3." The boundary conditions of diffraction are: arcs AB and 
A'B' that are the fronts of perturbations produced by small angle e, 
wall AOA', and curved shock wave BCB' resulting from diffraction away 
from the vertex. 


The problem is to determine the parameters of motion in region (3). 
Since angle € is small, quantities 


» Dp = a - 
Pu a 2 o=%—79, GC. VU) == 93—Ge (1.6) 





i) / 


2° 2 


will be of order e€, and the theory of perturbations can be used to find 
these quantities. Since there is no characteristic dimension of length, 
the problem can be considered self-similar, and then the equations of 
gas dynamics after linearization with respect to ¢€ in variables 


are written as 


Ov, a 0: Oy, = Oy, 


On plane &, n, the diffracted shock wave BCB' and the wall AOA' can be 
approximated by straight-line segments BB' and AA' (Fig. 2). 
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Fig. 2 


In this context, the value of §& on the wall is equal to zero, and on the 
curved shock the value of & is 


— bo-9" Ly (byl) 
Cy 





Arcs AB and A'B' in system &, n coincide with the arc of a unit circle. 
The solution of system (1.7) should satisfy the defined conditions on 


the boundary of the diffraction region. On perturbation fronts AB and 
A'B', the functions are continuous, and consequently on these arcs 


P=(), p=, g(u,v) =0 (1.8) 
According to (1.2), we have on the wall 


SPy Pat tateP) Pe gp 
Mgt Mgls 


implying that on the wall in region (3) 
” 2 aie 


Then, as implied by the second equation of (1.7), at € = 0 


OP _ 4, OF (1.9) 

0: 0", 
To find the boundary condition on & = Ky, we represent the perturbed 
shock wave front in the form — = k, + W(n), where p(n) is a function 


that defines the unknown shape of the shock wave. Writing the condi- 
tions of dynamic consistency on the curved part of shock BB' and 
linearizing them relative to the flow behind the plane reflected wave, 
we get the following system of equations: 
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where the projections on the normal and the tangent to the shock wave 
are written in the linear approximation in terms of projections on axes 
Ox, Oy, V,. + V' is the velocity of the curved shock wave, V' = 


cy(p ~ p'). 


The solution of system (1.10) takes the form 
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P = Qvy—12 ) E, u = G(y—y7v )/E, v=—q,» —49,(e +3) 
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Eliminating function p from (1.11), we get 
u=zA,P, ny on %2=Ry (1.12) 
Ov 07, 
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where A, = G/Q, Bo q,E/Q. 


Now by using equations (1.7) we can eliminate u and v from equations 
(1.12), giving the condition on the shock for P alone 
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In addition to condition (1.13), a condition is taken according to which 
the change in v along the shock from the center C to vertex B is 


(* p) 
—- dar, = — — dr=y: i.14 
. Ov J t OF ( ) 


- f — y= 9 (1.15) 


The Chaplygin-Buseman transformation is used to solve the formulated 
boundary value problem with respect to the solution of equation (1.15) 
with boundary conditions (1.8), (1.9), (1.13), (1.14): 
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converting (1.5) to a Laplace equation in variables p, 8B. Then confor- 
mal mapping is used [2]: 


ft -+ . ) 
sy eat —_— ): em TE § ( x nt, y= Xt ly) (1.17) 











Fig. 3 


Function (1.17) transforms diffraction region AA'B'B to a rectangle 
(Fig. 3) 


pal (1.18) 


Thus, the problem is to determine the function P that satisfies the 
Laplace equation AP = O in variables x,, y, within rectangle (12.18), 
and the conditions: 


gr 
on AB andA'B’ we have = —— == (1.19) 


OX, 
OP/Ox . 
AA’: ———- =%cot y (1.20) 
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2 oP 9 R \} , 2\105 
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In virtue of equation (1.20), there is a difference between the given 
problem and the problem of [2] in which @ = 0, while the fairly complex 


function w(z,) = se - ge is determined by trial and error. 


§2. Solution of Formulated Boundary Value Problem by Method of Separa- 
tion of Variables 


We will seek the solution of the Laplace equation AP = O in the form [3] 


Pix, y=) YO) 


In this case we get 
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X= Acxplzx,) + Bexp(—«x,), )=Csinxy, -Oeoszy, 
Satisfying condition (1.19), we get 
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ke) 
Condition (1.20) can be written as 
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and condition (1.21) can be written as 
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Multiplying (2.2) and (2.3) by cos ly, (1 = 0, 1, 2,...) and integrating 
from -t to 1, we get a system of equations for determining A, and B,. 
It can be shown by simultaneous solution of these systems of equations 
that A, and B, with even k are equal to zero, while A, and B, for odd k 
will be determined from the following system of equations: 


A(1—5)=B8,(144), (1+ 5)A,—(1 9) 8, =3] (1 - 9)A,- (1 6) 83] 
i “ a a Bok 
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R=2(n+1 \(Ao~2 Tok (Ana rex p((n ~1).)+B,.1exp(— (n+ 1)-)] 
0 
N =(n+1)[(1+49)Ang:.—(1—9) Bais] 
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From (2.4), all A, and B, are determined in terms of A, (nm > 2). Condi- 
tion (1.14) is used to find A,. In what follows, A, = A,/A,, Be = 
B./A,. The convergence of series (2.1) is proved as in [3]. 


§3. Determination of Gas (Fluid) Flow Parameters in Diffraction Region 


We write expression (2.1) for perturbed pressure in variables £&, n. 
From formulas (1.15) and (1.17) we have 
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Components u and v can be determined from the equations of system (1.7). 
To find coefficient A,, we find component v on & = O from the third 
equation of system (1.7) 


U= | a had (—1)"7 (Qn + WC Aan et +B oat) 
- Oy, aan') 
~ | 
“| arcsinyn-- SW —s'n(2k arc sin 4) (3.3) 
ra R 


The quantity v at € = 0, n = 0 is equal to the value of q.., and hence 


ALS (VD (204 1 Anne + Bars) = 29. 


n=!) ~ 


where q., is determined from formula (1.4). Hence we have 
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15 











A, nae 
> ( 1)" t(Qn+1)a(Aansi b Bonat) 


ne 


Numerical calculations are given to clarify the qualitative picture of 
reflection of a shock wave from the wall. From formula (3.2), the 
pressure on the wall is computed at Mach numbers M = V,/cyo = 1.2, 1.6, 2 
and @ = O, 0.1, 0.3. Calculations for the distribution of dimensionless 
pressure P = -P/eU, on the wall show that the pattern that occurs for a 
rigid wall is qualitatively retained, and moreover there is a decline in 
pressure P with decreasing 6. The distribution of total pressure P. on 
the wall is also calculated for e« = 0.1, pressure P. decreasing with 
increasing pressure @. The table also shows that P and P., on the wall 
in the diffraction region increase with an increase in the Mach number 
at fixed @. 


1) 

: 5 

05 | 
1 9- 

1 | 8-0 








Q= Of  -95 
“05, 8-03 ) 
-{ 
-{ | 7 
M=1,2 M=46 
Fig. 4 


Fig. 4 shows graphs of the distribution of pressure P on the wall for 
two cases. 
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[Text] In the linear theory of bending of plates under a transverse 
load, it is known that forces T,, (i, k = 1, 2) do not arise in the 
middle plane of the plate, and hence it is assumed that the transverse 
load acting on the plate does not affect frequencies of free oscilla- 
tions. However, in bending of very thin (flexible) plates, considerable 
tensile and shear forces arise in the middle plane of the plate for 
flexures comparable with the thickness of the plate, and these forces 
may significantly alter the frequencies of free oscillations of the 
plate. 


This paper attempts to give a qualitative picture of the influence of 
transverse loading on the frequencies of free oscillations of a plate. 


1. Let a flexible rectangular plate of dimensions a, b, h be hinged 
around the edge and subjected to the action of a uniform transverse load 
of intensity g. 


The equations of bending of a flexible plate referenced to a rectangular 
cartesian coordinate system Oxyz are [1] 





a... 2a) 32 ~, 3 
| y@, Pv a (<2 ) =z () 
ii dx? ay" OXdy » 

(1.1) 


m 2 ) y P 2 


0x? oy? "Ox dv Ox dy Oy? ax? 





where wo, % are functions of flexure and forces of the initial static 
state, D = EH*/12(1 - v*), E is Young's modulus, v is Poisson's ratio. 


17 








the forces in the middle plane are defined as tollows in terms of func- 
tion %,: 


72 2 2rt) 
T" = @ D, T° _— Avy T" Oo [ 0 


i ; : : (12) 
ay" - Ox? | dx OY 





Small transverse oscillations of the plate about the bent equilibrium 
state (1.1) are described by the equations [2] 








Azp + Pw Orie PW Oe 9 O°, ew oO 
Eh Oy? Ax? Adx* dy? 0x Oy AXaY 
Pw , Pwo Pw ., Ow 
ah DMT! — —2TY, -T), — - (1.3) 
Ot 0x* OX Oy “ Oy 
Av, a°cb 5 Cui, a, LT LAS | 
eisai — —— == ( 
dy 0x* JC Ov AX OV ax* ay 


where w, @ are the functions of flexure and forces of the perturbed 
(oscillatory) state, p is density of the plate material, and t is time. 


Equations (1.3) account for both the forces TY, and flexures w, of the 
initial state. 


Equations 
Dy=Noqy, Slur 7 ye stu ys C'y == Foxy SIE WV Stn yy (1.4) 
= Darp SIN XN SIN MAY COS Waal, Wfinr SLi aN SIN AV COS OmaAl (1.5) 
where iA,, = mt/a, u,, = nt/b identically satisfy the boundary conditions 


for hinged support of a plate 


/ =), w=0,. «=f, M,, as) for Xz), Xa 
(1.6) 


Tog 20), uN, we—--N, M,,=9 for v7), y= 0 
As shown in [1], expression (1.4) within the limits of flexures of 
initial state w, of the order of two or three plate thicknesses gives a 


rather good description of the initial bent state. 


For (1.4), (1.5), using the Bubnov-Galerkin method, systems of equations 
(1.1), (1.3) yield 


6 
Ky fou al n= = / (1.7) 


(Won ny = 


, oo 
’ ma dn f>) (1.8) 
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Here we have introduced the notation w,,, to denote the frequencies of 
free oscillations of the plate about the bent equilibrium position, and 








ie 
“oma=|/ m9 (40 rit) (1.9) 
to denote the frequencies of free oscillations of the loaded plate 
. , , Eh ; — : 
A mnt im * -)?, Ann = — | —“, - oe 2 | i2mn 
: 3 F207 - ey 
(1.40) 
fag ee | Dine 
| att? (4? —1) (47-1) | }/m fo yttn 





As can be seen from (1.8), the frequencies of free oscillations of a 
plate loaded by pressure q increase as compared with a plate free of 
pressure, and clearly the influence of the initial state can be disre- 
garded for flexures of the initial state that are much less than the 
thickness. 


2. By way of illustration, consider the example of a plate of square 
planform at v = 0.3. In this case, equation (1.7) is written as 





Prose Y.=O0 (2 ") 
where f, = fo51,/h, q, = 9a*%q/8n°h*E are dimensionless parameters of 
deflection and loading. 

TABLE 1 
Sn Sn nie i a onl Mens Reali 
qe | | 2 | 3 4 WW | 1S | 20 
j | | ee 
7 | a sn es 
24 


| 0.37 | 0.67 | 0.99 ey ee ee ee 2 





Table 1 gives values of f, for several values of q,. 


From (1.8) for frequencies of free oscillations of a loaded square plate 
at v = 0.3 we have 


A 


‘ 
1s) = ee 
mir 2 
a 


‘ “1m? n=)? { (1 San 2 f* 7 (- ~) 


where we have introduced the notation 


mint(in®? —n?—O0,9)_ 
om" (mt —0),25)(n?—0,25) 
On the basis of (2.2) for various values of q, we determine the normal- 
ized frequencies of free oscillations of a loaded plate: 


” 


a” 


Oy? (2.3) 











for m = n= 1 for which w,,-— min. 


m,n 


TABLE 2 
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Clearly, accounting for forces and flexures of the initial state has a 
fundamental effect on the frequencies of free oscillations of the loaded 
plate for flexures comparable with thickness. 


The numerical results obtained in section 2 are general for arbitrary 
plates of square planform. However, it can easily be shown that these 
results are competent only tor sufficiently thin plates that allow 
bending comparable with the thickness of the plate within the limits of 
elasticity or up to brittle fracture. From the condition 


9? Gs, Tee $ 35) +9519 (2.4) 
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“99 





1 ab, Ez (Ze uly lad, Ez Ow 

h Ox" [-.* \ dy? Ox? / : i Oxdy 1+ Oxdy 

we can get the following estimate for the admissible relative thickness 
h/a: 


| adm 14,17. —9,85972) 


Here f, for prescribed q is determined from equation (2.1). 


Thus, for prescribed a, E, O,am, there are real plates for which ac- 
counting for the initial state is significant in determining the fre- 
quencies of free oscillations. For example, at O,4,/E = 0.005 and f, = 
2.1, a = 200 cm, we get h $ 2.53 cm. Let us note that f, = 2.1 corre- 
sponds to a normalized load of q, = 15. 


In conclusion, we note that questions of the influence of an external 
statically applied load on the frequencies of free oscillations have 
been rather thoroughly studied. And in the case of shells, the initial 
state caused by external loads may be treated as zero-torque [3] or 
Significantly torqued, as well as nonlinear [2, 4]. On the other hand, 
in the case of plates, only geometric nonlinearity enables accounting 
for the effect of the initial state on frequencies of small oscillations 








about the bent equilibrium state. In this sense, we can derive rela- 
tions for studying the frequencies of free oscillations of plates of 
rectangular planform and for circular plates loaded by transverse pres- 
sure from the results of [2, 4] obtained for a cylindrical panel [2] and 
for corrugated shells of revolution [4]. 


Let us also note that we can derive an approximate analytical expression 
on the basis of [2] for the frequency of free oscillations of a loaded 
plate, and [4] can be a basis for exact numerical analysis of frequen- 
cies of a circular pressure-loaded plate. 
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[Text] Consider the problem of designing an orthotropic circular plate 
of radius R and thickness h under constant normal load q. 


The plate is referenced to a cylindrical coordinate system 0, r, 6, z2 so 
that coordinate plane z = O coincides with the middle plane of the 
plate. The material of the plate has three planes of elastic symmetry 
that are perpendicular at each point to the corresponding coordinate 
lines, the pole of anisotropy coinciding with the center of the plate. 
The plate of greatest stiffness is determined by varying the structure 
of the material and the location of the support contour. 


It is assumed that the plate is made up of 2n layers alternately rein- 
forced in the radial and circular directions. In this arrangement, the 
layers with odd numbering (1, 3,..., 2n-1) are made up of m, elementary 
layers that are reinforced in the radial direction, while the layers 
with even numbering (2, 4,..., 2n) are made up of m, elementary layers 
reinforced in the circular direction. From the foregoing we see that 
the total thickness is h = nmé,, where 6, is the thickness of a mono- 
layer, and m = m, + m,. The parameter &€ = m,/m defines the relative 
thickness of the radially reinforced layers of the stacked plate. It is 
assumed that elementary layers reinforced in both the radial and circu- 
lar directions have the same content of reinforcing material per unit of 
volume, so that the elastic characteristics of these layers in orthogo- 
nal directions are the same [1]. 


Let the plate be supported around line r = R,, where O < R, S$ R. The 
equation of equilibrium of an orthotropic plate relative to bending W(r) 
is represented as [2] 























where i = 1, 2 for regions 0 Sr 


7 
ry 
WA 


SR and R, R respectively. 


Boundary conditions are written as 


) 


(3) 





are the bending moment and the transverse force in cross’ sections 
r = const respectively, and D,,, Deg and D,. are the stiffnesses of the 
plate of given structure [1] 






(A = Bats 
D, lee npinecinaae 


I _ l | a f. 
Here B,,, Bee and B,,. are the elastic characteristics of a monolayer of 
orthotropic composite material. 


In deriving these expressions, terms of order 1/n and 1/n”* were ignored, 
where n is the total number of monolayers reinforced in the radial and 
annular directions. 


The solution of boundary value problem (1) and (2) with allowance for 
(3) takes the form 
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where the notation 
» _ Dw Biy—1): : | 
pee OT rR, HR,IR 
D,, (1—Bia)i = By, 
DA B,. 
=— =, 
D,, ( \ Ki ) 1 By, 


has been introduced, in which B,. B,.e/B,, and Bee = Bee /B,-,-- 


Our problem is to find the optimum value of parameters € and 


maximize the stiffness of the plate, i.e. 
f — iil iia \ . ) 
Woot witht ’ W _ 
with constraints 
Oil, Uscsisel, Os (il), Seayscl (/=2) 


and for assigned R, h, B,,, Bee, Bre and q. 


B that 


(S) 


By way of example, consider a plate made up of monolayers of orthotropic 


boroplastic reinforced in the radial and annular directions. 


An algorithm was developed for numerical solution of problem (7) with 


constraints (8), and the corresponding program was written. 




















| — 
: i W’ 
| 
0,00 0,693 100 06452 
0,05 0,691 1 -UU 0.5463 
0,10 0.689 1 QU U 4834 
0:15 0.688 0,00 0 +4325 
(+20 0,687 0.00 0.3896 
0425 0,685 1 00 03688 
0>30 0.684 1-00 0.3437 
0,35 0683 1,QO 0 +3235 
0.40 +682 10 0.3072 
0.45 0,681 1,00 0.294 
(655 0.660 | 0.00 "2674 
0.60 U.679 | 0.00 VU +2550 
0,65 0.678 | 1.00 0. ‘486 
Or? 0.677 1.00 0.2442 
0.75 0,676 | 1 ,Q0 0,410 
0.80 0,676 | 0.00 0.2332 
0,85 0,675 0,00 0.2262 
1,90 0,674 0,00 0.2211 
0,95 0,673 | 0,00 0,2195 
100 0,672 0.00 0.2262 
Values of Bupe = B,, the value of y = r/R where max W,(y, 8, &) is 
reached, and the corresponding min max W,(y, 8, &) are given in the 


i] Y 
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table for several values of structure parameter € characterizing the 
anisotropy of the plate. Here 


Was W 648,49 / log 4 
is is the dimensionless flexure of the plate. 


As we see from the table, the optimum value of parameter 8 that charac- 
terizes the position of the supporting edge for all values of structure 
parameter & varies monotonically in the narrow range 


0672. 3% 0,693 
The optimum plate design is realized for 
50,99, 40,673 and here \t opt= 2195 


this design containing 5% of monolayers reinforced in annular direc- 
tions. The worst design is realized for &§ = 0 (plate made up of 
monolayers reinforced only in the annular direction), and here 


Inax minimax Wes) 645z 


Thus, structural optimization leads to considerable improvement of the 
design of a plate in the sense of stiffness, 


maka aan a eR? Gt antnominmas ou 


i.e. roughly by a factor of 3. 
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[Text] The paper examines the axisymmetric problem of elastic wave 
propagation in a radially inhomogeneous elastic ball with the assumption 
that radial and tangential stresses on its surface are known. Standing 
oscillations of a radially inhomogeneous elastic ball are also studied. 


1. General formulas. Let the Lamé coefficients A, yu and density p of 


an elastic ball be assigned in spherical coordinates (r, 8, @») by func- 
tions 


ren (J a wy(—) * = (>) (1.1) 


lo Fay u/ 
where Ag = (Y - 2)Ho; Hor Po are the Lamé coefficients and density of 
the medium at r = rg; Y = (Ao + 2Uo)/Uo- 


The wave field of perturbations in an inhomogeneous elastic ball r < ry 
at any time t > O is determined by some equation of elasticity 


(A+ 2p)grad div u—s curl curl «@—grad 1. div a + 





scwnldh + Sated'n prada (1.2) 
+ gradp  curla + Z2grad eg er yt 
under conditions 
i<- «0 (1.3) 
Ot |r ~0 


Srey, 9 t)=Bayf(8, E)i slur 8. C)=ZoFal%s O); Zre(Tor se =f #1) O14) 
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Assuming Rayleigh's hypothesis, according to which the relation between 
components of the stress tensor and displacement vector U(u,, Us, u,) of 
an inhomogeneous elastic medium is the same as for a homogeneous medium, 
we have in the axisymmetric case 


Ou | Ou, ~ Ole us _ 
i= On — ’ rr ( : ~ = (1.5) 
Or r of Or r 
1 ou ll, Ou u 
ot 1) (-) — Qn ( - — ip it ( — *) ( l 0) 
[ Oh r Or A 
7 | / Ou u 
3 = 4 ' ( cot 1) Th, tt l ans noee +f ( Le/ ) 
where r r r of 
| a 
A= a sll (rae oe (sin [us ) (1.8) 
r* or rsin4 ot 


The displacement vector U of the radially inhomogeneous medium is repre- 
sented in terms of generalized potentials as follows: 


\ 


f prX\@= 


ulr, 9, (1 grad Dr, 6, N+(= our curl (Cr, 6, t)e,|+curl[D(r,8, te] 


/ P 
{ ; 


(1.9) 


i.e. the components of the displacement vector ti are determined from the 
formulas 








ab a a): 
u, ad at ES (eins ) (1.10) 
Ok R sin mer 05 | 
l aw* ine aD, 
enh exe ment -!. ° ee (1.11) 
w=r| Rap * ako | 
Pes lid, (1.12) 
” R A 
where @. is the basis vector with respect to the r-axis, R = r/Yro, %% = 


o,/r2, i= 1, 2, 3. 


Then, representation (1.9) of the displacement vector permits separation 
of the vector equation of motion (1.2) into the following independent 
scalar equations [1]: 


4 ad: 1 ; 1 @icb: 
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where 
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2. Wave propagation in an inhomogeneous elastic ball when its surface 
is acted on by radial time-dependent pressures. In this case, gen- 


eralized potential %* depends on R and t, o* = @* = 0. Wave propagation 
in the ball is described by the equation 





BJ)" +e ‘p" , 1 ap) e 
me 5 sansoniatinange oti. s . - 0 4) 
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under conditions 

0) Seas) 

Saas 2. ‘ 

— 7 p' | . —3 | | a “0 f(<) (2.3) 
0R* 1 tr 


Equation (2.1) and condition (2.3) in the Laplace transform with respect 
to T 


re . 2. 
4.1(R, p)= | DCR, Dexp(—pz)d= (2.4) 
are as follows: 
ao, 2d a : si 7 (2.5) 
dR® ° (-—2)R dk ~  O) R? 
: _. (42 5 
yo aD, 7-2 aD, ir; 9 
v wy For ) ane >_—_—- { ) (2.6) 
-R ape 2)R aR |. me p 


The solution of Euler's equation (2.5) under condition (2.6), assuming 
boundedness in the center of the ball, takes the form 


tof (p)R*-° 








* =— 9, 
(RR, ~P) “"g(7—C)(7— a) (2.7) 
where 
~ - —_ 4 / > — 
( 2 —9)* d=c— 7 ; 1= f/ ad c? — 
The radical branch a is fixed by the condition arg a = O for p > 0. 


According to (1.10)-(1.12) and (1.5)-(1.8), The components of the dis- 
placement vector and stress tensor take the form 
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where 


aTIVA| PR’ lor! 


ata ve, (2—d) 
s+(R, p) sof Dyk 
9 /) 
til Ry p) 3, (R. p)=- Kr (2 , \R 
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(2.8) 


(2.9) 


(2.10) 


We get the inverse Fourier transform of images (2.8)-(2.10) under condi- 


tion f(p) 


1, i.e. f(t) = 6(t), where 6 is the Dirac function. 


In the 


general case for any f(p), we find the inversion of expressions (2.8)- 
(2.10) from a convolution formula. 


From formulas [ 2} 


we get 
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i | 
where H is the Heaviside function, J, is a Bessel function. 


Waves of displacements and stresses propagate in a radially inhomo- 
geneous elastic ball when its surface is acted on by pressures that 
depend only on time. The front of these waves at any t > O is defined 
by the equation 


ee (2.15) 


Equation (2.15) defines a sphere with radius R = exp(-v yt). As a result 
of rapid decay of the radius of the wave front in time, the displacement 
and stress waves have an appreciable effect near the surface of the 
ball. These waves do not reach the center of the ball, and consequently 
no reflected waves originate from the center. 


In formulas (2.12)-(2.14), the terms associated with the Heaviside 
function describe waves that do not originate in a homogeneous medium. 
It propagates with a wake, and its occurrence can be attributed to 
dispersion of transmitted waves in the inhomogeneous medium. Near front 
(2.15), the so-called wave with a wake has a finite discontinuity im- 
plied by formulas (2.12)-(2.14) if it be considered that the limit 
z °J,(bz) as z > O is finite. 


Let us consider steady-state waves in a radially inhomogeneous ball. 
For boundary condition f(t) = o,e'”* we get 





Veo inRl] va (2.17) 


The phase velocity is determined from the condition 2 = const: 
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The wave number is 





. ‘)) l a 7 4 . ) | 9) 
A ) ) S—lClU _ —— So aa aes ies bas 
( U BR | 7 J 1 


For the group velocity we get the expression 





oa oe | Pg (2.20) 
As we see from formulas (2.18) and (2.20), the phase and group veloci- 
ties depend on frequency w, coordinate R, and the quantity y. 


Standing waves appear in a radially inhomogeneous ball at frequencies 
lower than the limit 


w=) (e+) (2.21) 


3. Resonant frequencies of a radially inhomogeneous elastic ball in the 
case of standing waves. Consider the solution of system of equations 
(1.13)-(1.15), assuming that the generalized potentials take the form 


PR, 4, =P, (R)P(cos Der, sel, 2.3 (3.1) 


where P,, cos 8 is a Legendre function, and radial and tangential stress- 
es on the surface of the ball are equal to zero. 


Substituting (3.1) in (1.13)-(1.15) and (1.4), we have system of equa- 
tions 





F de, } Ww Ff, __ ) 
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3 4 4 5 72 _ 
Poy 8 ean ty. “0 rae Y (3-3) 
dk? (7 -2)R dR Lt _~ 

ae, dz, n(azly) +r — roe (3.4) 

dk* (,-2)RdR 1 ~* 
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where 








tp =3Pg Fr ~n(nt lk “| Pawns) (3.8) 
dR 
+: . | 
|i p uk aP,(cos?) (3.9) 
u= 0 k df a4 
ry Paivosh) (3.10) 
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The solution of Euler equations (3.2)-(3.4) under condition of bounded- 
ness in the center of the ball is as follows: 


o(R)=AnR™"*; O(RI=BR" 's 45(R)=CrR” (3.11) 


where A,,, B,,, C, are arbitrary constants, 


———— —_———————— 


. —_—cooC ’ ” 
/ 7 iu f b- “2s ‘, - Cf 
- —_—- / e = ome , n d 4 
i ) a —_—_— | 3n=) n | 1 
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Radial oscillations of an inhomogeneous elastic ball occur at n = O. 
Solution (3.11) at n = O with boundary condition (3.5)-(3.7) gives the 
equation of frequencies 


4 . 
ena dy a? . mrsen(n Vy) ch mtap es ats ty eh 


n . —-\) oan 


The roots of equation (3.12) that define the resonant frequencies of 
radial oscillations are as follows: 


i29 -« (>, —a) - 0 (3.12) 


Roots of frequency equation A,,(w) = 0 




















| 0.0016 | 1.2245 | 1.6942 252565 2,9270 3.6715 
255 2.3665 | 28647 | 4.0385 | 4.4722 4:4722 | 4.4722 4.4722 
4.4722 | 4.4722 | 4.4729 | 5.4739 6.9535 | 8 +3606 
9.0030 | 1,3040 1.9568 | 26994 | 3.4642 
$4642 | 3.4642 | 3.4642 | 4.1854 54777 | | 
| | | 
0.0083 | 1.4912 | 2.2362 | 2.2352 | 2.2362 | 2.2362 
1U 24.362 | 2+2362 | 2.2362 | 2.4089 3.3412 4.2849 | 5.2334 
352849 | 2,6439 | | | 
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At n = O, substituting (3.11) in boundary condition (3.5)-(3.7). we get 
a homogeneous system of linear algebraic equations relative to A,, B,, 
C,,. The condition under which the homogeneous system has a nontrivial 
solution is as follows: 


S() =4(2, —O 0,0) ~~ —4]—(~--Qyn(ns p(s, - oN3,—-e—2) — 


Syn Wy) dq - 1)(a,--¢ -1)(4,-¢ - 3) =9 (3.14) 


The roots of equation (3.14) define the dimensionless resonant frequen- 
cies of a radially inhomogeneous elastic ball. They are given in the 
table for several values of n and y. 
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Introduction of Rigid Cylinder Into Cylinder Made of Material Conforming 
to Power Hardening Law 
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Sciences | 


[Text] In [1-3], researchers consider introduction of a cylindrical 
solid into an ideally plastic tube (isotropic and orthotropic cases). 


Our paper takes a look at introduction (internal and external) of a near 
cylindrical solid into a cylinder made of material with power-law hard- 
ening. 


1. Let a cylindrical tube made of material taken as incompressible and 
power-law hardenable fit tightly in an absolutely rigid cylindrical 
mold. The inside and outside radii of the tube are a and b respective- 
ly. A cylindrical tube of considerably harder material with varying 


outside radius R = R(z) = a + vu, exp (vz/b), where v and u, are pre- 
scribed positive constants (Fig. 1), is coaxially forced into this tube 
(internal introduction). The material of the second tube is taken as 


nondeformable. 


The problem is solved in the cylindrical coordinate system r, 6, z fixed 


to the rigid tube so that plane z = O passes through the entrance end 
cross section, and the positive direction of the Oz axis is opposite to 
the direction of motion along the axis of the tube. The end z = 1 of 


the tube is assumed to be free of external forces. Then the components 
of displacement u, w in directions r, z do not depend on 6, and v= 0. 


The equations of equilibrium, relations between the components of 
strains, displacements and stresses, and the law of hardening take the 
form 
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We seek the components of displacements in the form 


4 1 
ry | m P i. - e\p( 2) .* +) 


r 


where f = f(r), w = (rf)', A is a constant, and the prime denotes the 


derivative with respect to r. 


Then | 
: | _), ——evpt *) (1.5) 


eXP(7l), 2 = — eNpl 








where 


From relations (1.2) and (1.5) we get 


- t . . , y on, _ 
sy = aye (J — = Jexp(mi->). 3 ont aeliie ak —) (Pp (m=) 
fr, 


trz2 =0,54,07 exp(mz) (1.6) 


where w = w(r) = e§™ '?7". 


From (1.1) and (1.2) we get 
ty) ae, . 
1, = H(z)--- Rk, exp(me-z) jo - T+ (—) )dr (1.7) 


From relations (1.1), (1.6), (1.7) we find the unknown function H(z), 
and for function f we get the fourth order differential equation 


H(2z)=12,C, exp (m-) 


/ ‘ ) _ \ / wir! 
wl) —mitw( mT + 2(— — 2mit(w(— ~r") oes : () 1.8 
ory —mita(nt + 2(Z))—ania(= =) =(F) <9 am 
where C, is a constant. 


Boundary conditions take the form 


tre(, Z)=4,3,-(, 2), t,2(0, =)=43,(8, =) (1.9) 


(1.10) 


Moreover, the condition 


u(6, z)=0. ula, Z)=u, exp(-s b) 


where uw, and uw. are coefficients of friction. 


| 22(7, Yyrdr=0) (i.U1) 


must be satisfied on end z = l. 


We introduce dimensionless variables 


r z 
uo= — 10, “s—, y=. =m, 9 — 
b b b h D 
u(r, 2)==bu*(o, .), wir, =) = bw tcl, SPS Ro, 5) (1.12) 
3p(7, =)==Ry3"(0, 5), se(r, =) Ssh (y, 3), 22(7, 2) = Ay2(% 5) 


fir)=b7*(2), R(z)—bR*(-), Cy=oc, 
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Fourth order differential equation 


(1.8) can be replaced by a system of 
four second order differential equations. 
it takes the form 


In dimensionless quantities, 











df* pty d=* 4.72 2* 
conan = = as ee * = = 
ue Jn ad. a J 
dk < ff . 
— = i fi 9° } (1.13) 
do ° j 
df* l ji* ,of ** + 4%" 
— = —(3F* | —_— - : —_ — (| — mM )+*5( 28," a 3)7*o*/ 
uo Y* 4 4.2 
5 —3m)2f 8T* (3 —m) 22287" 4 42m —1) 405" FF 
— 987 . +/7l*" Sad ft ** arg *?) 
where 
3 ° f*) 
.* . —*)y, *® i _ 7 a —_ 
‘ 
r s 272 / e2 
»% (e*) . — ), e* ¢ #—|— : - _ 
40- + 
(1.14) 
Fe = fy f- (T*m*) — 2Qartg*.* T tus 
Yeah Qa ee) Em) 
Using relations (1.6), (1.7), (1.12), (1.14), we can find the dimension- 
less components of stresses 
:, (i : 3,\ byexp ii ) { ( * ‘ : Vie ‘yl ) 
. = (C* 7“ > 3. ) ), = iy { /?i } \ >) 
where 
iW) = 
(1.16) 
() | ' 
—— */ 3 ly * ! . | --.7)* 
; i l | | \ } } ) } 
The value of C* in (1.16) is obtained from (1.11). 


Boundary conditions and dimensionless quantities take the form 
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P*(4, )w* (9) == T*(ljo*( l= (CC, (i)), 


2. At a small value of parameter v we set 





where 
f=fotfit?h— 2B, re, we ee Te 
Fes Thole 0 gt oy Foe FytoFyr OF gt 
e= Li pr : : ye ee Rot ue (2.2) 
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an as, ) ’ z a] ‘ ‘ —y - 
= Cyt Cy tC ...5 we, 


Considering relations (2.2) at v = 0, (1.13) and (1.17) give the system 
of equations 


dfy _ MtunSy, Go _ Sho= 0%, WF yg 


2." do 


7) a ds 
“ | | (2.3) 
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and boundary conditions 


rif 1)=0; Ly(%) = Uy, Tyo) g(%) — 24Ciy 
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3fy— 49 
Ty(Lag (1) = 5 (20 | oa) 
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Solving the first three equations of (2.3) and satisfying the first two 
boundary conditions of (2.4), we find 


; | ; as Ne 
So= pie " :) 9, = —C, (3+ <) Py= 20, (.) 
where C. = PoUo/(1 - pg) and C, is an arbitrary constant. 


Substituting (2.6) in the last equation of system (2.3) and making 
necessary simplifications, we get a nonlinear differential equation for 
To 





whose integral takes the form 


clue) Pa (orep)™ e 
a a COOP 


where C., is an arbitrary constant. 





Constants C, and C,, are determined from the last two conditions of (2.4) 
which, after substitution of f, and 9, take the form 


+ , + a 42 
C+ Cite — toh f (. aot do 


) a ’ 
l ~ Or 4 
e/ af 





(2.9) 
I 
~ . ‘y ( 7 “4 “Pr 
C, C- | (3 | AOE de 
] - I "4 le 
At m = 1, we find from (2.8) 
; C , Cc ? 
[= (2.10) 
and conditions (2.9) take the form 
Ce Cy iCute t tur C,+ Cy =~ 8Cay (2.11) 
At m = QO, we get 
| ps Cy (2.12) 
ry = 2.) 3 V4 (C5 Cy" 
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The problem has been numerically studied at py, = 0.3, pw, = 0.2, a = 10, 
b = 12.5, u, = 20, v = 0.025 (uy = 1.6, Po = 0.8) for different values 
of m: m = 0.1, 0.2,..., 0.9. Values of P*/{, have been calculated, 
where P* = P/(k,b”), P is the force of pressing 


| ee 4m... \ 
' =” " ’ - oom ‘, ) "4 
vs ac | K*) | te * Pras ray’ ( “~" * ( My ‘) )} ul v) yl itu! 


0.6F : 





Lt i : i 
) 02 04 0608 1m. 
Fig. 2 





Fig. 2 shows how P*/1 depends on the quantity m. 


3. In the case of external introduction (Fig. 3) R* = 1 - VWuUyg exp (Vv) 
and the boundary conditions take the form 


Tehihs > )='1,3-(K, =). tr2(0, - = 1,7,(0, z) ( 1) 
u(d, z)= lyeXpP(-20), la, 2) —U (3.2) 


Formulas (1.1)-(1.8), (1.11)-(1.16), (2.1)-(2.3) hold in this case, 
while the remaining formuias change 
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ELECTRICAL TRANSIENT PROCESSES IN THE VAPOR-GAS SHELL UNDER CONDITIONS OF 
ELECTROLYTIC ANODIC HEATING WITH RECTANGULAR VOLTAGE PULSES 
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[Article by S. I. Galanin, V. I. Ganchar, and E. G. Dmitriyev, Institute of 
Applied Physics, Moldavian SSR Academy of Sciences] 


[Text] Anodic electrolytic heating (AEH) is accompanied by film boiling of 
the electrolyte solution on the surface of the anode. The two-phase system 
which results is electrically conductive, and the existence of a vapor film 
(shell) around the anode results from the joulian heat release of the current 
flowing through it. One of the methods of obtaining information about the 
nature of the current carriers in the gaseous phase is analysis of the 
dependence of the electrical conductivity of the vapor-gas shell on the 
applied voltage. However, such dependence cannot be obtaiied directly from 
the volt-ampere characteristic of the stationary process of AEH, for in 
stationary conditions the change in electrical resistance of the shell under 
change in the voltage results primarily from change in the thickness of the 
shell [1]. The required dependence of the shell conductivity on the applied 
voltage may be obtained by pulsed measurement methods. In the most simple 
experimental layout for a pulse technique, it is sufficient to measure the 
time constant of decay of the discharge current of a capacitor of given 
capacitance through the shell. Such a simple method is only useful when the 
shell represents a purely active resistance. In reality, when the electric 
field is turned on in pulses, the phase boundaries of the regions near the 
electrodes are charged and dual electric layers are created on them, which may 
be approximated as capacitors having a certain overall capacitance C,. In 
order to clarify the influence of this capacitance on the measurements of 
electrical conductance of the shell, it is necessary to know its magnitude. 
This is the task of the present work. 


In order to measure the magnitude of the overall surface capacitance C,, 
rectangular voltage pulses were used. A steel anode was heated in a 

stationary AEH process up to T 2 800°C, and then the voltage was switched off. 
A single rectangular voltage pulse was applied to the thus-prepared current- 
free shell in a time much shorter than its persistence time. The results of 
the action of this pulse on the electrolytic cell with and without an anode 
vapor shell present in it were qualitatively different (see Fig. 1): the 

electrolytic cells does not distort the original shape of the imposed pulse; a 
distortion is only observed when the current-free anode shell is present in 
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the system. This experimental fact suggests that the presence of the current- 
free anode shell in the electrolytic cell brings about an additional 
capacitance in the system, sufficient to distort significantly the shape of 
the imposed pulse. 

















( 
a 
500 742). 








Fig. 1. Comparison of current and voltage 
oscillograms obtained with a rectangular 
voltage pulse acting on an electrolytic 
cell in the presence (1) and absence (2) 
of a current-free vapor shell in it. 

Key: a. us. 


It is easy to support this conclusion with calculation. For this, we repre- 
sent the electrolytic cell with current-free anodic shell present in it as an 
equivalent substitution circuit (Fig. 2). The current and voltage measured by 
the oscillograph may be obtained from the system of equations: 


UL. =/y(R, T R, tT Ros) +] pRn, 
] I 

— fledi=h-R, ” 
Cc 0 

I, =In+Ic, 1p(0)=0, 








where U, is the potential of the source of rectangular pulses; I, is the 
current flowing through the electrolytic cell, I, and I, are the currents 
flowing through the active resistance of the phase interface surfaces R, and 
the overall surface capacitance C,, respectively; Ry, Rg, and Reg are the 
internal resistance of the source and the active resistances of the electro- 
lyte and shell, respectively. From (1), we obtain for the sought quantities 
U(t) and I,(t): 











vaya , Bly es Ramet yop] 2 I), 
(t)=e "= RR; | — R» CaRa ) (2) 

fo(ty = Cf ft exp _ =x) 
R,; 3 \ R> aRn (3) 
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mare R,=R,+Roo, R.=R,+ Root Ru, R3=R,+ Root Rar 


Calculation by (2) and (3) yields time functions U(t) and I,(t), similar to 
the experimental. The observed distortion in shape of the imposed pulse is 
achieved when C, > 2000 uF for values of R,, R,, and R, on the order of 
0.1-1 ohms and duration of the imposed pulse r = 500 us. 











| | 
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Fig. 2. Equivalent electrical substitution 
circuit of an electrolytic cell with 
current-free vapor shell present in it. 


The surface capacitance of the phase interface boundary can also be evaluated 
from other considerations. As an example, consider an electrolyte-vapor 
interface surface. Under the action of the imposed voltage, the phase 
boundary is charged until the electric field of the surplus charge accumulat- 
ing on the surface cancels out the external field. The characteristic 
dimension of the charge separation will be the Debye screening radius, which 
is connected to the equilibrium charge density within the electrolyte n by the 
familiar relation: 

lp =(eT/4 re2n)!/2, 


(4) 


Here, T is the temperature of the electrolyte near the surface, expressed in 
energy units; e is the charge of the electron; « is the dielectric 
permittivity of the electrolyte near the phase boundary. The surface layer of 
electrolyte, rp) thick, may be approximated as a flat capacitor, the capaci- 
tance of which is: 


Ca _ @&/T p, 


(3) 


where ¢, = 8.85-10"* F/m is the dielectric constant. As a result, we find that 
C, ~ /n, and as the electrolyte concentration varies from 13 to 38 percent the 
specific surface conductivity of the electrolyte-vapor phase boundary changes 
from 0.82 to 1.5 F/m*. In the conditions under which the measurements were 
taken, the anode surface was 2.7-10° m’; thus, the total capacitance of the 
phase surface should change from 2200 to 3900 uF as the electrolyte 
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concentration varies in the indicated range. Such dependence of the surface 
capacitance on the electrolyte concentration is experimentally confirmed: as 
the electrolyte concentration rises, the shape of the imposed voltage pulses 
is more heavily distorted (Fig. 3). 
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Fig. 3. Oscillograms of voltage at various concentra- 
tions of electrolyte NH,NO,: 1 - 32 percent, 2 - 
13 percent. Solid lines - calculation, points - 
measurement results. 
Key: 

a. V b. us. 


As for the nature of the influence of the surface capacitance of the phase 
boundaries on the pulse measurement results, it may be variable. For example, 
when low-voltage pulses are applied to the current-free anode shell and the 
potential U, is less than a certain critical value U,, there are no conduction 
currents through the vapor shell [2], and the surface capacitance C, in the 
equivalent substitution circuit is connected in series with the proper 
capacitance of the shell Cgg. Since Coé<< C,, the overall capacitance of the 
system is C ~ Cgg. Thus, the influence of the surface capacitance is 
eliminated, making it possible to measure the proper capacitance of the anode 
shell. Such measurements were performed in [3]; according to the data of this 
work, Cg9g= (150-300) pF. A different situation occurs when measuring the 
conductance of the anode shell: when a capacitor charged to potential U, > U, 
is connected to the shell, a discharge current flows through it; the 
capacitance of the phase interface surface is added to the capacitance being 
discharged. The capacitance of the phase interface surface is on the order of 
thousands of microfarads, while the electrolyte-vapor phase surface is able to 
discharge in a time on the order of hundreds of microseconds. Consequently, 
in order to eliminate the influence of the capacitance of the phase boundary 
on the measurements, the discharging capacitance must be chosen in the range 
of values much larger than the capacitance of the phase boundary, and the 
current decay curve must be registered in the millisecond time range in a 
period longer than the charge time of the phase interface surface. 
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Shamko, Electrohydraulic Design Development Office, Ukrainian SSR Academy of 
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[Text] The results of experimental studies of an electrical explosion in 
ferromagnetic suspensions [1, 2] show that, depending on the conduction and 
concentration of the solid phase, both the conditions of formation of the 
channel and the discharge regime vary. This should result in change in the 
gas dynamic characteristics of the underwater electrical explosion (UEE) and, 
accordingly, the effectiveness of transformation of electrical energy into 
mechanical [3, 4]. Therefore, the present work is devoted to developing the 
methods of calculation of the energy and gas dynamic characteristics of an 
electric explosion in a liquid containing ferromagnetic inclusions. 


In the first part of the work, the methods of physical modeling are used to 
investigate the electrical, energy, and kinematic characteristics of the UEE 
with initiation by ferromagnetic particles, localized in lengthwise bridges. 


The experiments were carried out in a tank (0.9 x 0.9 x 0.6 m’) with portholes 
of organic glass, filled with distilled water of electric conductance 

oo = 2.5010° &' x m' (the layout is similar to that of [5]). In order to 
reduce the cumulative effect of the electric explosion cavity (EEC), the 
negative electrode had a rectangular opening of 0.03 x 0.23 x 0.1 m (the 
external appearance of the electrode system with ferromagnetic bridge is shown 
in Fig. 1). The distance between electrodes in all cases was 1 = 40 mm, and 
the discharge occurred at depth of 0.3 m. 


The bridges were created from ferromagnetic powder SF-35 with initial electric 
conductivity o, = 10° Q'em' as follows. In a separate vessel, a 10 percent 
suspension of powder SF-35 was prepared from distilled water, and a capillary 
was used to supply this to the interelectrode gap (IEG) with a lengthwise 
magnetic field, created by a special coil [1]. When the suspension enters the 
region of the field, the solid particles line up along the field lines and 
form a bridge spanning the entire IEG (see Fig. 1). 
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Fig. 1. External appearance of electrode 
system and ferromagnetic bridge: 1 - 
negative electrode; 2 - positive 
electrode; 3 - ferromagnetic bridge. 


The electrotechnical parameters (charge voltage of the capacitor bank U,, its 
capacitance C, the inductance of the discharge circuit L) are presented in 
Table 1. For each particular set of parameters U,, C, L, at least 

13-15 discharges were performed, and the obtained results were averaged. 


Table 1. Electrotechnical parameters of the discharge circuit 








Key: a. Number of discharge 
b. V Cc. F d. J 


NOTE: L, 10° = 3.5 H. 


The studies were done by the method simultaneous registration of electric and 
space-time characteristics of the discharge channel and the EEC. The current 
was registered by a low-resistance coaxial shunt, the voltage drop across the 
discharge gap by a capacitive divider, the signals from which were sent to the 
S8-17 oscillograph. The inductive component of the voltage drop was compen- 
sated. The discharge channel was photographed by the high speed camera VFU-1, 
acting as a photoregister (mirror rotational speed 30,000 rpm), the cavity was 
photographed by the SFR-2M with an adapter reducing the mirror turning speed 
by two orders of magnitude. Highlighting during investigation of the channel 
stage was provided by the pulse lamp IFK-2000, and after the discharge the 
constant light source KG-110 was used. 
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The current and voltage oscillograms (similar to those presented in [1] were 
used to determine: in the channel formation stage, the breakdown delay time 
t,, the decrease in magnitude of the initial potential AU,, the energy losses 
AW (averaged values of these are given in Table 2 for the different discharge 
regimes), and in the active stage the energy generated in the channel W(t), 
the time for its generation r, and the electric power N(t). Deviations of the 
investigated quantities from the mean values usually did not exceed 30 percent 
for t,, 50 percent for AU,, 10 percent for r, and 15 percent for AW/W,. 


Table 2. Averaged time and energy characteristics 
of discharges 1-3 
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Key: 
a. Number of dis- b. us 
charge S. kV 


A characteristic of the formation of the channel when the discharge is 
initiated by powder SF-35 is the fact that, within 20-30 us after voltage is 
applied to the IEG, several luminous regions of size 1-2 mm are created at its 
center, the number of which increases in the lengthwise direction over the 
course of time (Fig. 2a, b). Occasionally, two luminous chains appear, 
spanning the entire gap at a certain instant. 


The mechanism of formation of the discharge channel (breakdown stage) when the 
discharge is initiated by conductive ferromagnetic inclusions has been 
discussed in [{l, 2, 6] and is explained by the redistribution of the electric 
field in the gap as a result of the presence of the bridge, resulting in 
concentration of current and joulian heat release in individual regions on the 
surface of the inhomogeneity. The latter, in the opinion of [1, 2, 6], may 
result in initiation of a discharge in these zones. It may be that a similar 
breakdown mechanism occurs in a su -ion based on the powder SF-35, since 
its electric conduction is greate’ that of water. However, other 
mechanisms may be at play (breakdv uiong the surface of the bridge, break- 
down of gas bubbles ever present in the liquid, etc.). 


It is difficult to distinguish in time and space the stages of formation of 
the discharge channel (prior to breakdown) and the channel proper, as is 
obvious from Fig. 2a. Toward this end, the discharge channel was photographed 
without lighting (see Fig. 2b) and the appearance of an intense diffuse light 
around the channel in the second frame was coordinated (on the analogy of a 
discharge without initiation) with the start of a sharp rise in current on the 
oscillograms. Thus, the electric and kinematic characteristics of the 
discharge when initiated by ferromagnetic bridge were synchronized. 
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characterizing the discharge regime and representing the portion of energy 
released during the first half-period of the current. Here, A is the spark 
constant, which for an electric explosion initiated by microconductor is 

A = 0.25010° V’es/m’. 








a,MmM a,m/c(a) 











0 20 4O 69 oMKU ) 








Fig. 3. Time curves of radius of channel 
a(t) and speed of its expansion 4(t) for 
the first discharge: 1 - experiment, 2 - 
calculation. 

Key: a. m/s b. ps 


In the present work, the parameter n was determined experimentally as the 
ratio between the energy W, released during the first half-period of the 
current and the stored energy W, = CU,°/2. It turned out that, for unchanged 
initial electrotechnical parameters of the channel and physico-mechanical 
properties of the working medium, the parameter changes from one discharge to 
another (for an UEE initiated by microconductor, the regime is uniquely 
characterized). The range of change in n for discharges 1-3 is shown in 
Table 2. While for conducting powders [1, 2, 6] this was due to change in 
concentration of the solid phase in the ferromagnetic suspension and, thus, 
the initial diameter of the bridge, for the nonconducting powder SF-35 this 
may be explained by an increase in length of the channel as compared to the 
length of the IEG (see the above expression for 7). Values of the equivalent 
channel lengths lg, corresponding to the UEE initiated by microconductor and 
calculated for different discharge regimes, are presented in Table 2. The 
above-given results of experimental investigations of the energy character- 
istics of an electric explosion initiated by SF-35 powder show that they can 
be determined in accordance with the technique [4] developed for an ordinary 
underwater spark discharge. However, in calculating the parameter n, instead 
of the length of the gap 1 it is necessary to use 21 < ly < 3l. 


52 











30> N,10°8m (a) 





#0 4 | 





4 Wi 


0 20 0 f,. 
4 "th) 


Fig. 4. Power curves for the first 
discharge: 1 - experiment, 2 - calcula- 
tion by [4]. 


Key: a. W b. us 











We investigated the influence of the opening in the negative electrode on the 
kinematics of the cavity. Since a poor quality of photograph of the EEC is 
obtained when the UEE is initiated by powder SF-35, Fig. 2c presents photo- 
graphs for the first discharge initiated by copper conductor of diameter 

O.1 mm. During the first pulsation, the cavity is dome shaped, and becomes 
slightly deformed in the closure stage. It has been established that the 
shape of the cavity, its maximum volume, and the period of the first pulsation 
are almost unaffected by the opening, which basically influences the speed of 
movement of the EEC boundaries (the cavity collapses more intensely if the 
opening is missing, see Fig. 5a, b) and the number of compression pulses 
generated by the cavity. Thus, when an opening is present, the EEC produces a 
single compression wave, and when it is absent (closed by a metal plate), it 
produces two waves. Analysis of the kinematic characteristics of the EEC 
(height of the dome dy(t), its speed dy(t) in the perpendicular direction (a) 
and the radius of the cavity a_(t), the speed ay (t) in the direction parallel 
to the negative electrode, Fig. 5) reveals that the functions more sensitive 
to presence of an opening are d,(t), d)(t). The kinematic characteristics of 
the EEC formed during UEE initiated by powder SF-35 for the third discharge 
are shown in Fig. 5c. 


In the concluding stages of the closure (see Fig. 2c), the cavity moves in the 
direction of the free surface and collapses to form a cumulation stream of 
liquid, 4 cm above the point of contact between the channel and the plate 
covering the opening. At the point of impact of the stream on the plate, a 
conical funnel is formed with radius of depression at the contact point 

1-2 mm. Because the regions of influence of the channel and the cavity 
(stream) are at a distance from each other, it was possible to separate their 
mechanical action on the plate. It was found that the bending of the plate 
(using copper and steel strips, 6 = 0.25-1 mm thick) at the point of action of 
the shaped stream of liquid was several millimeters larger than at the point 
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of contact of the channel. Thus, for a copper plate with 6 = 0.25 mm, the 
bending was respectively 10 and 7 mm. This result is consistent with the data 
of [7, 8]. 















































a,,mm 
/ 
/ 100 
50 
7 10 15 t,MC 
/ 3 
Ee 
Fig. 5. Kinematic characteristics of cavities: a, b - initiated 
by copper conductor, 0.1 mm diameter, first discharge; c - 
initiated by powder SF-35, third discharge; 1, 2 - experiment 
(primed numbers - radius), 3 - calculation; --- electrode without 


opening. 


Thus, the initiation of a discharge by nonconducting ferromagnetic particles 
satisfactorily stabilizes its energy, kinematic, and time characteristics. 
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